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Master Equations for Gravitational Perturbations 
of Static Lovelock Black Holes in Higher Dimensions 

Tomohiro Takahashi and Jiro Soda 
Department of Physics, Kyoto University, Kyoto, 606-8502, Japan 

In this paper, we derive the master equations for gravitational perturbations of vector 
and scalar type for static vacuum Lovelock black holes. Together with our previous work on 
the tensor type perturbation, we now provide the full set of master equations that governs all 
types of gravitational perturbations of static vacuum Lovelock black holes in any dimensions. 

t-H ■ 

§1. Introduction 

There is a long history in black hole perturbation theory since the seminal paper 
qj [ by Regge and WheelerP~® n j s we n known that there exist master equations for 

P-h ' gravitational perturbations of static black holes in 4-dimensionsP^ Surprisingly, 

the master equations have been obtained even for stationary black holesP'^ 1 It 
should be mentioned that the master equations for gravitational perturbations have 
i_— played an important role in gravitational wave physics. 

^ | Recently, higher dimensional black holes have attracted much attention. This is 

i | because higher dimensional black holes could be created at the LHC.^ In fact, in 

^Jy the context of the braneworld with large extra-dimensions, the predicted production 

rate of black holes is within reach of accelerators. Another reason is that higher 
\ dimensional black holes have been used to analyze strongly coupled finite temper- 

1 ature field theories through the AdS/CFT correspondence. Needless to say, the 

master equations of gravitational perturbations of black holes in higher dimensions 
^ | are crucial for the developments of these subjects. 

t— I ' The master equations for gravitational perturbations of higher dimensional static 

black holes have been obtained by Kodama and IshibashiP For stationary black 
holes, unfortunately, there exist only partial results. To investigate the stability 
of rotating black holes, a group theoretical method is developed.® The method is 
used to obtain the master equations for gravitational perturbations of squashed black 
.£h ! holes^ 1 "® and 5-dimensional rotating black holes with equal angular momenta.^ 

The master equations for a special class of rotating black holes in more than 5- 
dimensions are also studiedPHS* Still, it is an open issue if the master equations 
exist for general rotating black holes and black rings. 

We should notice that, in higher dimensions, Einstein theory is not a unique the- 
ory with the second-order differential equations. Indeed, the most general theory of 
gravity is Lovelock theory which is degenerated into Einstein theory in 4-dimensions. 
In fact, Lovelock theory is a natural extension of Einstein theory in that the Love- 
lock theory contains terms only up to the second order derivatives in the equations 
of motion. In Lovelock theory, it is known that there exist static spherical symmetric 
black hole solutions^ 1 and topological black hole solutions.^ 1 Hence, it is legiti- 
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mate to suppose black holes produced at the LHC are of this typeP^ 1 Thus, it is 
important to obtain master equations in order to study these Lovelock black holes. 

The following argument also stresses the importance of Lovelock theory. The 
point is that black holes are produced at the fundamental scale of higher dimensional 
theories. At such high energy, Einstein theory would be no longer valid. In fact, as is 
well known, string theory predicts Einstein theory only in the low energy limit P In 
string theory, there are higher curvature corrections in addition to Einstein-Hilbert 
termP3 Thus, it is natural to extend gravitational theory into those with higher 
power of curvature in higher dimensions. It is Lovelock theory that belongs to such 
class of theories 

In the case of second order Lovelock theory, the so-called Einstein-Gauss-Bonnet 
theory, the master equation for tensor perturbations has been obtained p3 The re- 
sult has been also extended to the scalar and vector perturbations P2J Although 
Einstein-Gauss-Bonnet theory is the most general theory in five and six dimensions, 
it is not so in more than six dimensions. For example, when we consider ten dimen- 
sional black holes, we need to incorporate the fourth order Lovelock term. Indeed, 
when we consider black holes at the LHC, it is important to consider these higher 
order Lovelock termsP^ Hence, in this paper, we derive the master equations for 
gravitational perturbations of black holes in any order Lovelock theory, namely, in 
any dimensions. We have already derived the master equation for tensor perturba- 
tions of static Lovelock black holes in any dimensions. 25 ' In this paper, we derive 
master equations for vector and scalar perturbations of static Lovelock black holes. 

The organization of this paper is as follows. In Section [21 we briefly review 
Lovelock theory and static Lovelock black hole solutions. In section [31 we express 
the linear Lovelock tensor in terms of the perturbed Riemann tensor. In section HI 
for completeness, we derive previous results for tensor perturbationsPS 1 In section 
[5J we derive the master equation for vector perturbations. In Section El we deduce 
the master equation for scalar perturbations. The final section [7] is devoted to the 
conclusion. 

§2. Static Lovelock Black Holes 

In this section, we review Lovelock theory and introduce static black hole solu- 
tions. 

In Ref. I21|) . the most general symmetric, divergence free rank (1,1) tensor is 
constructed out of a metric and its first and second derivatives. The corresponding 
Lagrangian can be constructed from m-th order Lovelock terms 



where R\rr pK is the Riemann tensor in ^-dimensions and 5^ 1 ° 1 '''^ m ° m is the general- 
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ized totally antisymmetric Kronecker delta defined by 

/ «5£ 6% ••• 6% \ 
6% 5% ■■■ C 



V S% 5% ••• J 
Then, Lovelock Lagrangian in D-dimensions is denned by 

k 

L = c m C m , (2"2) 

m=0 

where we defined the maximum order k = [(D— 1)/2] and c m are arbitrary constants. 
Here, [z] represents the maximum integer satisfying [z] < z. Hereafter, we set 
Co = —2A, c\ = a% = 1 and c m = a m /m (m > 2), for convenience. Taking variation 
of the Lagrangian with respect to the metric, we can derive the Lovelock equation 

k 

n — C, v — AR V — \ am yAiffi-A m g m p p\n\ t> Pm^m (n o\ 

u ~~ y y. ~ yu V Z^- 2( m+1 ) m AtpiK1 "' PmKm 1 " A " lff '» ' 



m=l 



As shown in Refs. [TT|) and I18p . there exist static exact black hole solutions of 
the Lovelock equations. Let us consider the following metric 

dr 2 



ds 2 = -f(r)dt 2 + — - - + r 2 j ij dx i dx j 

fin 



(2-4) 



where jij is the metric of n = D — 2-dimensional constant curvature space with a 
curvature k=1,0 or -1. Using this metric ansatz, we can calculate Riemann tensor 
components as 



r> tr 
Ktr 



ki 



f" tj rj f j 

' i Rti — Rri — 2r 

K f\ ( t* k c I c- Z c- & 

—2- Ibi bj -Si bj 



Substituting (|2-5p into f|2-3|) and defining a new variable ^(r) by 

fir) = k — r 2 ifj(r) , 

we obtain an algebraic equation 



m=2 



in 



a 
m 



2m-2 



II in-p)\r 

p=l 



+ V 



2A 



n(n + I) 



(2-5) 
(2-6) 

(2-7) 



In (|2-7|) . we used n = D — 2 and is a constant of integration which is related to 
the ADM mass of black holes as^ 1 

M = JgL — - , (2-8) 

r((n + l)/2)' 1 ' 

where we used a unit 16ttG = 1. 

In the following sections, we study gravitational perturbations around the gen- 
eral vacuum solutions (|2-6[) obtained by solving algebraic equation (|2-7|) . 
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§3. Linear Lovelock Tensor 

In this section, we present general formulas for the linear Lovelock tensor around 
the solution (|2-4p . The method introduced in this section is based on symmetry of 
the static Lovelock black holes. In detail, the point is that only Rt r tr , Rtfi , Rri r ^ 
and Rij kl have nonzero value for the Riemann tensor due to the metric ansatz (|2-4p . 

From now on, we use //, u, ■■■ for (t,r,x l ), x, y, ■■■ for (r, x l ), and i, j, ■■■ 
for (x l ), respectively. With this notation, we can show that 

C£r^ = {n-(m- m%:::£ o-i) 

by induction. This formula is useful for later calculations. It is also easy to see that 
the linear Lovelock tensor reads 



k 

2(m+l) '• v-!>n, ih„ 

xR Xiai ^ ■ ■ ■ Rx m ^ m .^ Rm -^Rx m .J mKm ■ (3-2) 



xnu _ _ \ ^ Um xv\icri-"\ m a m 
y M Z_> o(m+l) i 

m=l 



Let us explain how to calculate 5Q\ in detail and merely present final results 
for other components. First of all, we consider the totally antisymmetric Kronecker 
delta in 5Q\, namely, ^tptKi'-'-p^Km '• Because of the antisymmetry, t cannot show up 
twice in the indexes. Hence, this Kronecker delta can be rewritten as 

r*2/l2/2 tt -J/2m-l2/2m _ r2/l2/2-"2/2m-l2/2m /q q\ 

°txiX2— X2m-lX2m ~ ° x l x 2- ■ ■X2m-l x 2m ! l '°7 

where we used properties 5\ = 1 and JJ. = 0. Thus, we can rewrite 8Q\ as 



yiV2—y2m-\y2m 

X2---X2m-\X2m 



v r> X1X2 p X2 m -3I2m-2ID I2ra-13:2m 

^ ll yiV2 ll y2m-3y2m~2 U1 l y2m-l J/2m 

Here, when taking the summation of x p and y p , we have to consider three cases. One 
is that there is no r index, that is, all x p and y p are i p and j p , respectively. For the 
other two cases, we have r index. Assume x p = r, then, there must exist p such 
that y p i = r because the totally antisymmetric delta consists of Kronecker delta. 

One possibility is x%q-i = r or X2 q = r (1 < q < m — 1). For this case, p' must be 
p = 2q — 1 or 2q because of the background property ()2-5[) . The remaining possibility 
is X2m~i = r or x<im = r - For this case, we have to take p = 2m — 1 or 2m. The 
reason is as follows. If p 7^ 2m — 1 nor p ^ 2m, there must exist q (1 < q < 2m — 2) 
such that y 1 = r. For example, if we take y\ = r, x\ or X2 must be r because of the 
formula (|2-5p . In any case, ^I^'..'. must be zero because of the antisymmetry. To 
summarize, 5Q\ can be written as 



SGt 



/ j \ 2 m +i 

m 

rjl-"i2mp. . »1«2 Y Y P. • i2m-3«2m-2^D. . *2m-l«2n 

il-*2m A JU2 ^J2m-3J2m-2 UJX 32m-132m 
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vP . "1 P. . J 2*3 v V R • *2m-4*2m-3XR . . *2m-2«2m-l 

AJTy-jj -Ttj 2 J3 A A Jrt J2m-4J2m-3 W - rt J2m-2 J2m- 1 



+ 45; 



■31— 32m-l 
U-»2m-l -"-JU2 



«1*2 



V D. . *2m-3*2m-2^p . *2m-l»" 

A JX J2m-3J2m-2 U IX 32m-\r 



(3-4) 



Substituting the background quantities (|2-5p into (|3-4|) and using the formula 
we can proceed as 



Ml = E 



a 



2m+i 



m=l 

ora-1 f K ~ f 



m—l 



S3l—32m Wl f 1 2m-2rn. . «2m-li2m 

2 / °il-i2m°il J 2m _2 0ja J2m-lj2m 



+4(m - 1)2 



m-2 



2r , 



m-2 



v rJl"-J2m-l sh A* 2m - 3 AP. . *2m-2*2m-l 

A °n-i2m-l °jl U j2 m -3 OIl n™-2j2m-l 

/ f \ m—l 



E< 

m=l 



2m- 1 



-J fj (n-p) 



p=2 



m—l 



2m-2 



p=2 
2m-2 

n ( n - ^) 

p=i 



4t" / \ v A 

m—l 



b\bR jr l 



2{n - l)r n ~ 2 1 k 3 r"" 1 1 3 



(3-5) 



where we used the relation b^bR^ = 2bfb l jbRki tJ in the last equality. Here, T(r 
is defined by 



T{r) = r^d^Wty] = r n - 1 1 + ^ 



m=2 



2m-2 

<U„ \ \{{n-p)\^ m - 1 

P =l 



(3-6) 



Similarly, we can deduce other components of the linear Lovelock tensor. The 
results are as follows; 



2(n - l)r r 
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SGI 



T^5R it ir , 

r.n—1 l J 1 ' 



SG\ 
SG r r 
SGI 
SG] 



»n-l 



(n — l)r" * r 

t' 7 1 

— ; s oSiSiSR, Al 

2(n - l)r n ~ 2 1 k 3 

T' „„• ^ T 



n— 1 « J 1 ' 



(n l)r» ' " 

t' 



SlSRr/ k + —SRrt* 



(n — l)r n 2 
T 



(SRa tj + + 



(n - l)(n - 2)r™- 3 



r.n—1 



5Rtr tT + 



(n — l)r n 2 

t" 



SRtk tl + SR r k rl ) Si 



2(n-l)(n-2)r n -^ 6 P 6RkqlP 



(3-7) 



Thus, in order to derive the linear Lovelock equations, we need to know the perturbed 
Riemann tensor SR U u pX - 

Since Lovelock black holes have n-dimensional symmetric space, we can classify 
perturbations into tensor, vector and scalar type perturbations. In following sections, 
we treat these three type of perturbations separately. 



§4. Master equation for tensor perturbations 



In this section, we derive the master equation for tensor perturbations of static 
Lovelock black holes. 

Tensor perturbations around the solution (|2 4|) is characterized by 

Sg a b = , Sg ai = , Sgij = r 2 4>(t,r)hij(x l ) , (4-1) 

where a,b = (t, r) and (j>(t, r) represents the dynamical degrees of freedom. Here, 
tensor harmonics hij are defined by 

y^khj = -Jthj , V%3 = , fihij = , (4-2) 

where V* denotes a covariant derivative with respect to 7y and the eigenvalue is 
given by jt = £(£ + n — 1) — 2, {I = 2, 3, 4 ■ ■ ■ ) for k = 1 and positive real numbers 
for k = —1,0. Note that indexes i, j, ■ ■ ■ are raised or lowered by 7^. 

For tensor perturbations, from (|3-7p . it is clear that components other than SQf 
vanish and the terms proportional to 5j in 5Gj also disappear. Therefore, what we 
have to calculate are SRtfi , 8R r i r ^ and b™b~RirJ l . From the metric ansatz ()24p and 
(|4-ip . these components can be deduced as 
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+ 



4 r 



srsRiJ 1 



n- 2 f , 2k + 7 t 
' 2 r 2r 2 



(4-3) 



Then, substituting these results into (|3-7|) . we can calculate as follows; 
(n - l)r n - 2 <5ft j 



rT 



re - 2 

/ _/ 



fT i fT | fT 



+ 



(2k + 7f )r" 

2(n - 2)r 



V ■ (4-4) 



Separating the variables 4>(r,t) = x( r ) e wt ■> we can derive the master equation 
for the tensor perturbations from the linear Lovelock equation 5Q^ = as follows; 



{2K + lt )f T" 
(<n-2)r T 



tX = u X 



(4-5) 



Here, we should stress that we have assumed nothing for Lovelock coefficients and 
/(r). Hence, the master equation we derived is quite general. 

Furthermore, introducing a new function •F(r) = xi r ) r \/T' (r) and using tortoise 
coordinate r* which is defined as dr*/dr = l//(r), we can transform the master 
equation (j4-5|) into Schrodinger type equation 



d 2 V 



dr 



*2 



+ V t (r(r*))W = iv 2 & , 



where we have defined the effective potential 



Vtir) 



(2 K + lt )fd\nT 



+ 



1 



f-rif-rrVr 



(n — 2)r dr ' rVT TJ dr \ J dr 
Here, we have assumed T' > in order to avoid ghost instability.^ 



(4-6) 



(4-7) 



§5. Master equation for vector perturbations 



In this section, we consider vector perturbations of static Lovelock black holes 
and derive the master equation. 

We take the Regge- Wheeler gauge 





( o 





Vi 










Wi 




y sym 


sym 






where Vi and Wi satisfy transverse condition V l Vi = V l itfj = and "sym" means 
symmetric part of metric perturbations. 
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From the ansatz ()5-ip . it is clear that 5Q\ = bQ T r = bQ\ = 0. The components 
of SR^^ which is necessary for calculations of non-zero components of 5Q^ are as 
follows; 



SRtk* 3 



SRrk 3 



5Rt r tl 



5Rt/ j 



1 

2r 3 



(V; ( v 1 ' -w % - -v l j - S l k (v j ' — up — -v> 



f ■ f 
2r 3 2r 2 



5Rtr lr = ~^V l + ^-r 



1 



+ 



2r 4 L 



/ i / i 



1 i i' -i ^ , 

f - iir — u 



+ V 3 \i ) + — [Wi^ + w j u . . 



4^2 



It) 



(5-2) 



where we use \i as a covariant derivative with respect to 7y instead of V,. Then, 
substituting these results into (j3-7|) . we can get the linear Lovelock tensor as follows: 



2r n+2 Sgi 
2r n+2 5Gi 



2(n - \)r n bQl 



n — 1 

r' 
/ 



/ 



{fT 



where we used the relation 



{rf +2( K -f)}T=(n+l)n , 



(5-3) 
(54) 
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and its derivative with respect to r. Then, after expanding metric perturbations by 
vector harmonics Vi, which satisfies V^V^t^ = —j v Vi with 7^ = £(£ + n — 1) — 1 
(I > 2) for k = 1 and non-negative numbers for k = 0, —1 , we can get the linear 
Lovelock equation as 

- k) vT' - f [rT (v'-w- f «) }' = , 
'^-k\wT' -?f(v -w-^vX = , (5-5) 



n-1 

r 
/ 



T v + ( fT'w ) = . 



Note that only two equations among these equations are independent. In fact, we 
can get the third equation in (|5 5[) by combining a derivative of the first equation 
with respect to t and a derivative of the second equation multiplied by f(r) with 
respect to r. Here, we use the second and third equations in (|5-5|) . We have not 
considered the exceptional mode 7^ = (n — 1)k since its treatment is well known . 

©,[2711 

In this paper, we assume T' is always positive otherwise there exists ghost insta- 
bility for tensor perturbations. Under this assumption, we can eliminate v from the 
second equation in (|5 5[) using the third equation. The resultant master equation is 
given by 

rT - - ? ( -4 (WV + ( - *y «, = . (5.6) 



/ / [r 2 T /J V n - 1 

Furthermore, under the assumption of positivity of T', we can introduce a new 
variable 

/ . /rri' 

X = -VT w . 

r 

We also perform Fourier transformation \ = ^e~ lbji . Thus, using the new variable 
and a tortoise coordinate r*, we can rewrite Eq. (|5-6p as 

- + V v {r)V = lu 2 V (5-7) 

where we have defined the effective potential for vector perturbations 

V v (r) = rV¥fd r (fdr^=) + - «) ^ . (5-8) 

Note that we have not assumed anything for f(r) and Lovelock coefficients a m except 
for the positivity of T . In this sense, this Schrodinger type master equation is quite 
general. 

§6. Master equation for Scalar Perturbations 



In this section, we derive the master equation for scalar perturbations of static 
Lovelock black holes. 
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We use the Zerilli gauge: 





/ fH 


Hi 





hfjii/ — 


I sym 


H/f 







\ sym 


sym 





3-1) 



where "sym" represents symmetric part of metric perturbations. 

Following Ref . I23|) . we derive the master equation from 5Qj = 0(i ^ j), 8Q\ = 0, 
Sgf = 0, 5Gi = and 8Q r r = 0. Taking look at Eqs.fHETj), we see that we need SR ti tj , 
5R r i r:) , SRij kl , 5Ru rj 5R r ii k and 5Rt r tl • With the Zerilli gauge, these components 
can be calculated as follows: 



5Ru tj 



6Rri rj 



-Ljf,^' +\--h 1 + —k - f -K' + -La 1 + L H ) si , 

2r 2 ' I r 2/ 4 2r 2r J 1 ' 

-J-HJi + f !(/£?)' - f -K' - f - K ' - t K ") Sj , 



5fl 



1 

2^2 



+ 



—H-—k + -k' + — k i # , 

2r 4 2 2r i 1 ' 



SRt r tl 



1 

2^2 



2/ r f 



where |i means a covariant derivative with respect to Ty. Then, we can calculate 
the linear Lovelock equations. Expanding metric perturbations in terms of scalar 
harmonics Y which satisfy VfcV fc y = — j s Y with 7 S = l{l + n — 1) for k = 1 and 
positive numbers for k = 0, —1, we obtain the linear Lovelock equations 



sg{ = o(i^j)^ t\h -h) = tt 'k , 



5Q r t = - 7s #i + n 



r(if - H) + r 2 K 



K 



, 



(6-3) 
(6-4) 



= <^ {-7 S T - w(/T)'} # - nTrfH + r2 / (n« - 7 s )if 



+ 



+ nf{r z T) }K +nTr z fK =0 



(6-5) 
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ggi = & rT [-H-K 



+T 



f 1 1 

— (H + H) + h'--H- -Hi 
2p } f 



(6-6) 



sa 



— 44> — : — ( — ill — K \ + rT [UK — "f s )K + nT r / H — 



/ 



-nr(Tf) H + 7s T# - nTr/i? = 



(6-7) 



Now we derive the master equation for scalar perturbations from these equations. 
We do not consider exceptional gauge dependent modes j s = and 7 S = tik modes 
since treatment of these modes is well knownP'^ We also assume T' > outside 
the horizon again. First of all, we show that all perturbative variables H, H±, H and 
K can be expressed by a single master function <j) defined later. From Eq. (|6-3p . H 
can be expressed by H and K as 



H = H 



rT" 



-K 



(6- 



By inspecting Eq. ()6-4[) . we see that it is convenient to define a master function <p as 



Hx 



f 



+ K 



(6-9) 



Then, we can express H using eft and K by integrating (j6 4|) with respect to t. The 
result is given by 



H 



Is 

— f l 

nj 



+ rK 



A{r) 
2^7 



K 



where 



A(r) = -2nf + 2^ s + nrf 



(6-10) 



(6-11) 



Note that there may exist an arbitrary function of r as a constant of integration. 
However, this function can be absorbed into the definition of <f>. From Eqs. ([6-8p . 
(|6-9p and (|6-10p , it turns out that we need to express K by (f> in order to express all 
variables in terms of (j). Substituting (|6T0p into (|6-5p . we obtain such a formula 



K 



2 
A 



nrf4> + 7s + nrf 



T 



where we used a relation 

(AT)' = 2 ls T' + n{(rf - 2f)T}' = 2( 7s - uk)t' , 



(6-12) 



(6-13) 



which can be derived from (|5-4p . 
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Now, we are in a position to derive the master equation for the master variable 
i. From Eqs. (|6-6p and (|6-7p . we can make the following combination 



nfr x (l.h.s of dSHD) + (l.h.s of (Jg 

- ( 7s - n*;)rT'^ + (W - H) + - if ) .(6-14) 



Substituting (^8]) . (follj) . (jbMO]) and (^42]) into the above equation (pFT4"l) . we get the 
master equation for scalar perturbations 



/V + / 2 In 



.4 2 



r 2 /T' 



+ 



where we have defined 
/ 



nTr 2 



2 { -^-^)(jsrT + nrZfT' 



, 



n(r 2 /T ) 



Finally, we change the normalization of <j) as 



A 



rVT 



It is also convenient to move on to Fourier space as x = ^ e 



-iojt 



(6-15) 



(6-16) 



(6-17) 



Substituting 



these into the master equation (|6-15p and using a tortoise coordinate r*, we obtain 
Schrodinger type equation 



- c%. & + V a {r)V = uj 2 V , 
where we have defined the effective potential for scalar perturbations 
, (rNT)' 



(6-18) 



V s {r) = 2 ls f 



nNTr 2 



N'T' 



f N' 2 f T' 2 

-Ldr (fd r N) + 2/ 2 ^ - t dr{fdrT ) + 2/ 2 ^ + 2/ 2 — 



. (6-19) 



We assumed that T is always positive because tensor perturbation has ghost in- 
stability if this assumption is not fulfilled. However, except for this assumption, we 
have not imposed any conditions on the Lovelock coefficient a m and f(r). Therefore, 
the master equation we have derived is again quite general. 



§7. Conclusion 



We have succeeded in deriving master equations for gravitational perturbations 
of static Lovelock black holes. The results can be regarded as a generalization of 
master equations in Einstein theory derived by Kodama and Ishibashi. Of course, it 
is possible to extend our analysis to charged black holes and other black holes PS)J2S) 
Moreover, our result would serve a starting point for studying stationary black holes 
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in Lovelock theory. It is also interesting to consider Euclidean version of our results 
in the light of black hole thermodynamics.® 

There are many applications of our master equations. In an accompanying 
paper, using the master equations, we show that static Lovelock black holes with 
small masses are unstable in the asymptotically flat cases PD The application to 
asymptotically AdS cases is also interesting from the point of view of the AdS/CFT 
correspondence,^"^ [ n particular, in relation to stability of holographic supercon- 
ductors 123 J29 
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